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Abstract In a classical covering problem, we are given a set of requests 
that we need to satisfy (fully or partially), by buying a subset of items 
at minimum cost. For example, in the fc-MST problem we want to find 
the cheapest tree spanning at least k nodes of an edge-weighted graph. 
Here nodes and edges represent requests and items, respectively. 
In this paper, we initiate the study of a new family of multi-layer cov- 
ering problems. Each such problem consists of a collection of h distinct 
instances of a standard covering problem {layers), with the constraint 
that all layers share the same set of requests. We identify two main sub- 
families of these problems: 

• in a union multi-layer problem, a request is satisfied if it is satisfied 
in at least one layer; 

• in an intersection multi-layer problem, a request is satisfied if it is 
satisfied in all layers. 

To see some natural applications, consider both generalizations of k- 
MST. Union fc-MST can model a problem where we are asked to connect 
at least k users to either one of two communication networks, e.g., a 
wireless and a wired network. On the other hand, intersection fc-MST 
can formalize the problem of providing both electricity and water to at 
least k users. 

We present a number of hardness and approximation results for union 
and intersection versions of several standard optimization problems: MST, 
Steiner tree, set cover, facility location, TSP, and their partial covering 
variants. 



1 Introduction 

In the fundamental Minimum Spanning Tree problem (MST), the goal is 
to compute the cheapest tree which spans all the n nodes of a given edge- 
weighted graph G = {V,E). To handle the subtleties of real-life applications, 
several natural generalisations and variants of the problem have been considered. 
For example, in the Steiner Tree problem we need to connect with a tree only 
a given subset W oi k terminal nodes. In the fc-MST problem instead, the 



goal is to connect at least k (arbitrary) nodes. One common feature of these 
generalizations is that we need to design a single network. However, this is often 
not the case in the applications. For example, suppose we want to provide at 
least k out of n users with both electricity and water. In this case, we cannot 
design the water and electricity infrastructures independently: our decisions on 
which users to reach have to be synchronized. 

Consider now another classic problem, the Travelling Salesman problem 
(TSP): here we are given a complete weighted graph, and the goal is to com- 
pute the minimum-length tour traversing all the nodes. Again, several natural 
generalizations and variants of the problem have been considered in the litera- 
ture. Still, all of them deal only with the case where there is a single network. 
However, there are natural applications which do not fit in this framework. For 
example, suppose you want to visit a set of places (bank, post office, etc.), and 
you can use your bike and your car. Of course, you cannot just reach a place by 
bike, and then suddenly switch to your car (that you left at home). Your trip 
must consist of a tour by bike and another tour by car, which together touch all 
the places that you need to visit. 

The above examples show the need for a new framework, which is able to 
capture coordinated decision-making over multiple optimization problems. 

Our results. In this paper we initiate the study of multi-layer covering problems. 
These problems are characterized by a set of h instances of a standard covering 
problem (layers), sharing a common set of n requests. The goal is satisfying, 
possibly partially, the requests by buying items in each layer at minimum total 
cost. We identify two main families of such problems: 

• Intersection problems. Here, as in the water-electricity example, a request 
is satisfied if it is satisfied in all the layers. 

• Union problems. Here, as in the car-bike example, a request is satisfied if 
it is satisfied in at least one layer. 

We provide hardness and approximation results for the union and intersection 
versions of several classical covering problems: MST, Steiner Tree, (Non- 
metric and Metric) Facility Location, TSP, and Set Cover. (Formal 
definitions are given at the end of this section). We focus on the partial covering 
variant of these problems, i.e. fc-MST, A;-Steiner Tree, etc.: here we need to 
satisfy a target number k of the n requests. This allows us to handle a wider 
spectrum of interesting problems. In fact, for intersection problems, if fc = n 
it is suSicient to compute an independent solution for each layer. On the other 
hand, some of the union problems above are interesting also for the case k = n. 
However, the results that we achieve for that case are qualitatively the same as 
for k < n. 

For Intersection versions of A;-MST, A;-Steiner Tree, A;-TSP, fc-SET 
Cover, A;-Metric Facility Location, and fc-NoNMETRic Facility Loca- 
tion, we show that: 

• Even for two layers, a polylogarithmic approximation for these problems 
would imply a polylogarithmic approximation for A;-Densest Subgraph. 



We recall that the best approximation for the latter problem is 0{ni^^) 
[9] and finding a polylogarithmic approximation is a major open problem. 
Indeed, many researchers believe that a polylogarithmic approximation does 
not exist, and exploit this assumption in their hardness reductions (see, e.g., 

m)- . ^ ^ 

• On the positive side, we give 0{k^ ^' '')-approximation algorithm^ for these 
problems. 

Note that, in the single-layer case, the above problems can be approximated 
within a constant or logarithmic factor. Hence, our results show that the com- 
plexity of natural intersection problems changes drastically from one to two 
layers. 

For Union versions of fc-MST, fc-SxEiNER Tree, fc-TSP and fc-MEXRic 
Facility Location we show that: 

• The problems are ^^(log A:)-hard to approximate for an unbounded number h 
of layers. Furthermore, there is a greedy 0(log fc)-approximation algorithm. 
For the first three problems this only holds for the rooted version — the 
unrooted case is inapproximable. 

• There is an LP-based algorithmic framework which provides 0(/i)-approximate 
solutions. Furthermore, the natural LPs involved have f2{h) integrality gap. 

We remark that Union /c-Set Cover and Union fc-NoNMETRic Facility 
Location can be solved by collapsing all layers into one, and hence they are 
less interesting with respect to the goals of this paper. 

Related Work. To the best of our knowledge, and somewhat surprisingly, ap- 
proximation algorithms for union and intersection problems seem to not have 
been studied in the literature, with the notable exception of Matroid Inter- 
section. However, differently from our problems, Matroid Intersection is 
solvable in polynomial time |12) . 

The term "multi-layer" has been used before in the literature, but with a 
meaning different from ours. Most often it refers to problems related to VLSI 
design, where we are given several planar layers on which the circuit has to be 
built [IS]. It also sometimes refers to multi-layer models of communication net- 
works that are composed of different physical and logical layers of communication 
devices [35] . 

The idea of introducing multiple cost functions into one optimization problem 
is the main theme of muUi- objective optimization. Standard and multi-criteria 
approximation algorithms have been developed for the multi-objective version 
of several classical problems, such as Shortest Path [27I33I34I43] . Spanning 
Tree [20122134137] . Matching [6I7I20I34] etc. (for a survey, see [13]). One could 
view these problems as having several layers with different costs. However, this 
setting is very different from our approach. In fact, solutions in different layers 
of multi-objective optimization problems have to be exactly the same, and the 
goal is to satisfy some constraints on each objective. 



* The O notation suppresses polylogarithmic factors. 



Partial covering problems (also known as problems with outliers), are well- 
studied in the literature: e.g., fc-MST |3l4l5ll0ll7ll5155] . fc-TSP [HT], /c-Metric 
Facility Location [11'29;, and fc-SET Cover |30I40| . Their generalization on 
multiple layers is significantly harder, as our results show. Note that our Union 
/c-Steiner Tree problem generalizes all of the following problems: fc-STEiNER 
Tree (and hence fc-MST), Prize-Collecting Steiner Tree (see the proof 
of Theorem O, and /c-Set Cover (see the proof of Theorem [7]). 

Rent-or-huy [14115125141] and buy-at-bulk [19l21l24l25li^ problems can be 
seen as multi-layer problems where edge weights in different layers differ by a 
multiplicative factor. In contrast, weights of different layers are unrelated in 
our framework. Union fc-TSP has some points in common with multi-depot 
versions of TSP [32l36j : also in that case multiple tours are computed; however, 
their weights are measured w.r.t. a unique weight function. 

Recently Krishnaswamy et al. j31) considered a matroid median problem, 
where a set of open centers must form an independent set from a matroid. This 
can be viewed as a generalisation of a Union problem, however in the matroid 
median problem all the centers are in the same metric space. This setting is 
less general than ours as it does not allow for modelling a setting with several 
completely unrelated metric spaces. 

Preliminaries. In covering problems we are given a set U oi n requests, and a 
set S of items, with costs w : S ^ ^>o- The goal is to satisfy all requests by 
selecting a subset of items at minimum cost. We already defined MST, Steiner 
Tree, and TSP. Here, nodes and edges represent requests and items (with 
costs w : £' — > M>o), respectively. In the Set Cover problem, requests are the 
elements of a universe and items S are subsets 5*1, ... , 5,„ ofW. Any Si satisfies 
all the V G Si. Nonmetric Facility Location is a generahzation of Set 
Cover, where we are given a set T of facilities, with opening costs a : J- ^ R>0) 
and a set C of clients, with connection costs w : C x ^ M>o. The goal is to 
compute a subset A of open facilities such that X^/ex ~^ X^cgC ^("-j ^) 
minimized. Here w{c,A) := min/g^ ti;(c, /). We also say that c is connected to 
(or served by) A{c) := SLTgimnf^j{w{c, f). If connection costs satisfy triangle 
inequality, the problem is called Metric Facility Location. 

We can naturally define partial covering versions for the above problems: 
/c-MST, fc-STEiNER Tree, fc-TSP, fc-NoNMETRic Facility Location, and 
fc-METRic Facility LocatioinH. 

It is straightforward to define union and intersection versions of the above 
problems. In the rest of this paper, the number of layers is denoted by h, and 
variables associated to layer i have an apex i (e.g., w^, o\ etc.), whereas OPT 
denotes the optimum solution, and opt its cost. By N we denote the total number 
of requests and items (in all layers) . 

^ In the literature fc-NoNMETRic Facility Location often means that we are allowed 
to open at most k facilities, while here we mean that we need to connect at least 
k clients. Similarly for fc-METRic Facility Location. Sometimes /c-Set Cover 
indicates a Set Cover instance where the largest cardinality of a set is k, while our 
problem is sometimes called Partial Set Cover. 



Figure 1 Approximation algorithm for 2-layer Intersection fc-SET Cover. 
For a e {1, 2}, a is the other value in {1, 2} 
1: procedure SCl{k,U, ,w^) 



2: K ^ 0, ^ ^, A' ^ 9 

3: repeat 
4: for a=l to 2 do 

5: for all X e do 

6: for 6 — 1 to min(fc - IJfl, do 

7: Solve one-layer Intersection fe-SET Cover problem on layer a 

8: with universe X \ K and target b. 

9: Let {a',b',X') be the loop iterators which provide a solution {K',A') 

10: minimizing the ratio of cost C' to number b' of covered elements. 

11: K ^KU K', A"' ^ A"' U {X'}, A""' ^ A""' U A' 

12: until \K\ = k 

13: return {K,A^,A'^) 



By standard reductions, a p-approximation for the fc-MST problem im- 
plies a 2/9-approximation for /c-Steiner Tree and /c-TSP. Moreover, a p- 
approximation for fc-TSP gives a 2p-approximation for fc-MST. Essentially, the 
same reductions extend to the union and intersection versions of these prob- 
lems. For this reason, in the rest of this paper we will consider the union and 
intersection version of fc-MST only. 

2 Intersection Problems 

2.1 Intersection A;-Set Cover 

In this section we present our approximation algorithm for Intersection fc-SET 
Cover. We recall that in this problem we are given h collections , S"^ , . . . , 
of subsets of a given universe U, where : 5* — )■ K>o is the cost of subsets in the 
ith collection. The goal is covering at least fc elements in all layers simultaneously, 
at minimum total cost. 

The basic idea behind our algorithm is as follows. We consider any set X in 
any layer, and any number j < fc of elements in X. We solve recursively, on the 
remaining layers, the intersection problem induced by X with target j. The base 
of the induction is obtained by solving a one-layer Intersection fc-SET Cover 
problem, using the greedy algorithm which provides a (1 -t- lnfc)-approximation 
[30]. We choose the set X and the cardinality j for which we obtain the best 
ratio of cost to number of covered elements. Next, we include covered elements 
in the solution under construction, and the problem is reduced consequently. 

In order to highlight the main ideas of our approach, we focus on the special 
case h = 2, and we neglect polylogarithmic factors in the analysis. It is easy, just 
more technical, to extend the same approach to h > 2 and to refine (slightly) 
the approximation factor (See Appendix IA.3P . 



Theorem 1. There is a 0{^/k)- approximation algorithm for Intersection k- 
Set Cover on two layers. 

Proof. Consider the algorithm in Figure [TJ Its running time is polynomial, since 
SCI procedure calls the one-layer greedy algorithm O(Nk^) times. 

Let {0^,0^) C 5^ X 52 be the optimal solution, and let Kq Q {^seO^S) (1 
{Us(zc>2S) be any set of k elements in the intersection. For each element x £ Kq 
and layer i = 1,2, let us fix a set (x) S that covers x. We prove that 
at each iteration of the main loop C'/b' = d{opt/^k- \K\). This implies that 
the total cost of the constructed solution is bounded by Y^\Zo 0{opt/y/k — i) = 
opt ■ diVk). 

Let K := y/k — \K\. We consider two cases, depending on whether there exists 
a set X in the optimal solution that covers at least k elements of Kq \ K . 

Case 1. Assume that there exists 1 < a < 2 and X e O", such that for 
at least k elements x of Kq \ K we have (x) = X. Let us focus on the 
moment when our algorithm considers taking the set X. Obviously we have 
K < k — \K\, therefore our algorithm considers covering b := k elements of 
X. As the optimal solution does it, it may be done with cost opt, so the call 
to the one layer algorithm returns a solution with cost 0{opt). Hence we have 
C'/b' = d{opt/^k- \K\). 

Case 2. For each 1 < a < 2 and every X G O", at most k elements of Ko\K 
satisfy O" {x) = X. For each x£Ko\K, let w{x) := w^{0^ (x)) + w^{0^ {x)) 
be the sum of the costs of sets covering x in the optimal solution. We have 

2 2 

xeKo\K a=l XeC xeKo\K: 0'^{x)=X a=lXeO'' 

Thus there exists xo G Kq \ K such that w{xo) < k ■ opt/\Ka \ K\. If we take 
any a and consider the iteration with X = (xq) and 6=1, the algorithm 
computes a set of minimum cost Co < w{xo) covering a;o. We can conclude that 

□ 

The proof of the following theorem is in Appendix lA.3l due to space limits. 

Theorem 2. There exists a (4fc^~^/'' \o^^^ (k))- approximation algorithm for In- 
tersection fc-NoNMETRic Facility Location (hence for Intersection k- 
Set CoVERj running in iVoe*) time. 

2.2 Intersection fc-MST 

In this section we present a simple approximation algorithm for Intersection 
/c-MST. Recall that here we are given a graph G = {V, E) on n nodes, and h 
edge- weight functions w-^, . . . ,w^. By taking the metric closures of we may 



assume that G is complete. The goal is computing a tree for each layer such 
that J2t w'{T') is minimized and | fl^ V{T')\ > k. 

The algorithm is very simple: We consider a new metric w defined as a sum 
wie) :— J2i for each e e _E, and compute a 2-approximate solution of the 

resulting (one-layer) fc-MST problem using the algorithm in |18j . 

In Appendix I A. 1 1 we prove the following theorem. 

Theorem 3. The Intersection fc-MST algorithm above is Idk^^^^'^ -approximate. 

The analysis of the approximation ratio of the above algorithm given in 
Theorem ini is tight up to a factor 0{h) (see Appendix IA.1[) . 

2.3 Approximation Hardness 

This section is devoted to the approximation hardness of Intersection fc-MST, 
Intersection fc-SET Cover (hence also of Intersection fc-NoNMETRic Fa- 
cility Location) and Intersection fc-METRic Facility Location. We use 
reductions from the fc-DENSEST Subgraph problem: find the induced subgraph 
on fc nodes with the largest possible number of edges. The fact that partial 
coverage problems can be as hard as fc-DENSEST Subgraph is already known. 
Hajiaghayi and Jain [53] use fc-DENSEST Subgraph to show that a partial cov- 
erage version of the Steiner Forest problem has no polylogarithmic approxi- 
mation. In particular they introduce the Minimum £-Edge Coverage problem 
where one is to find the minimum number of vertices in a graph, whose induced 
subgraph has at least £ edges. Moreover Hajiaghayi and Jain show a relation be- 
tween approximation ratios for fc-DENSEST Subgraph and Minimum ^-Edge 
Coverage. In order to simplify our reductions we extend the result on Mini- 
mum ^-Edge Coverage to bipartite graphs and prove the following theorems 
in Appendix lA. 21 

Theorem 4. If there exists an f {n)- approximation algorithm for unweighted I'N- 
TERSECTiON fc-SET CovER On two layers or for Intersection fc-METRic Fa- 
cility Location on two layers, then there exists a I6{f {2m))^ -approximation 
algorithm for fc-DENSEST Subgraph. 

Theorem 5. // there exists an f (n)- approximation algorithm for INTERSEC- 
TION fc-MST on two layers, then there exists a 16{f {2n-i-2m+2))^ -approximation 
algorithm for fc-DENSEST SUBGRAPH. 

Theorems |4] and [5] suggest that the existence of a polylogarithmic approxi- 
mation algorithm for the considered problems is rather unlikely (or at least very 
hard to achieve). 

3 Union Problems 

In this section we present our results for Union fc-MST and Union fc-METRiC 
Facility Location. For Union fc-MST, we first consider the rooted case in 
Sections IST] and [3?2] and then the unrooted one in Section |3^ The Union MST 
problem is a variant of Union fc-MST with k = n. 



3.1 Approximation Hardness 



Theorem 6. Rooted Union /c-MST and Union /c-Metric Facility Loca- 
tion are APX-hard for any h>l. UNION MST is APX-hard for any h>2. 

Proof. The first claim trivially follows from the APX-hardness |18I23) of the con- 
sidered problems for ft, = 1, by adding dummy layers with infinite edge weights. 

For the second claim, we consider a reduction from the APX-hard [S] Prize- 
COLLECTING Steiner Tree problem: given an undirected graph G = (V^E), 
edge weights w : E ^ K>o, a root node r € V, and node prizes p : V ^ R>Oi 
find a tree T 3 r which minimizes X]eGT^(^) + Sd^tP(''^)- create a first 
layer, with edge weights = w. Then we construct a second layer, where we set 
w^{{r,v}) = p(y) for any v £V. All the other layers, if any, are dummy layers 
defined as above. This reduction is approximation preserving. □ 

For an unbounded number of layers, our problems become much harder. 

Theorem 7. For an arbitrary number of layers, rooted UNION fc-MST and 
Union /c-Metric Facility Location are not approximable better than i7(log k) 
unless P — NP , even when k ~ n. 

Proof. We prove the claim for rooted Union fc-MST, by giving a reduction 
from cardinality Set Cover: given a universe U of n' elements, and a collection 
iS = {5'i, . . . , Sm'} of m' subsets ofU, find a minimum cardinality subset ACS 
which spans hi. This problem is I?(logn')-hard to approximate [39]. We create 
one node per element of plus two extra nodes r and s. We create one layer i 
for each set Si (i.e., h = m'). In layer i we let w*({r, s}) = 1 and ui'({s, v}) = 
for each v € Si. We also let r"^ :— r for each i, and assume k ~ n ~ n' + 2. Note 
that any solution to the rooted Union fc-MST instance of cost a can be turned 
into a solution to the Set Cover instance of the same cost, and vice versa. 

To prove the claim for Union fc-METRic Facility Location, we use the 
same reduction as above, where the edge {r, s} is replaced by a single node r, 
which is a facility of opening cost 1. □ 

In Appendix IB . II we give a greedy 0(log fc)-approximation algorithm. 

3.2 An LP-Based Approximation for rooted Union fc-MST 

In this section we present an LP-based 0(ft,)-approximation algorithm for rooted 
Union fc-MST. Essentially the same approach works also for Union fc-METRiC 
Facility Location (see Appendix IB.2p . This is an improvement over the 
0(log fc)-approximation given by the greedy algorithm for the relevant case of 
bounded h. 

For notational convenience, we assume that the roots R :— Ui{r*} are not 
counted into the target number fc of connected nodes. In other terms, we replace 
fc by fc— . We make the same assumption also in the case of one layer. Consider 



the following LP relaxation for /c-Steiner Tree {W 9 r is the set of terminals) 
denoted by LP^sriw, W, V, r, k): 



min J2, 
s.t. E, 




\f{v,S) : S 'ZV - {r},v e S nW; 



> 0, 1 > Z„ > 0, 



Vw e wye e E. 



Here, variable Xe indicates whether edge e is included in the solution, whereas 
variable Zy indicates whether terminal v is connected. Moreover 6{S) denotes the 
set of edges with exactly one endpoint in S. Observe that LPkMSriw, V, r, k) := 
LPksriw, V, V, r, k) is an LP relaxation for fc-MST. We need the following lem- 
mas. 

Lemma 8. \17l Let {w, V, r, k) he an instance o//c-MST, Wmax ■= max^gy w{r, v), 
and opt' be the optimal solution to LPkMST{w,V,r,k) . There is a polynomial- 
time algorithm apx-kmst which computes a solution to the instance of cost at 

most 2opt' + Wmax ■ 

Lemma 9. 116128}/ Let G — (yu{w}, E) he a directed graph, with edge capacities 
a : E R>o such that J2ees+iu) ^i^) = J2ees-(u) "(e) Z"'' all u ^ VU {v}. 
Then there is a pair of edges {u,v) and {v,z), such that the following capacity 
reservation j3 supports the same flow as a between any pair of nodes in V : for 
Aa := min{a(u, v),a{v, z)}, set l3{u, v) = a(u, v) — Aa, l3{v, z) = a{v, z) — Aa, 
P{u, z) — a{u, z) + Aa, and /3{e) — oi{e) for the remaining edges e. 

Corollary 10. Given a feasible solution {x,z) to LPkST{w,W,V,r,k), there is 
a feasible solution {x',z') to LPkMST{w,W,r,k) such that X^e '^(e)^e — ' 

Proof. Variables Xe can be interpreted as a capacity reservation which supports a 
fractional flow of value z^ from each v G W to the root. Let us replace each edge 
with two oppositely directed edges, and assign to each such edge the same weight 
and capacity as the original edge. This way, we obtain a capacity reservation a 
which costs twice the original capacity reservation, and satisfies the condition of 
Lemma O We consider any non-terminal node v ^ r with some incident edge 
of positive capacity, and apply Lemma IHl to it. Due to triangle inequality, the 
cost of the capacity reservation does not increase. We iterate the process on the 
resulting capacity reservation. Within a polynomial number of steps, we obtain 
a capacity reservation /3 which: (1) supports the same flow from each terminal 
to the root r as a, (2) has value on edges incident to non-terminal nodes 
(besides r), and (3) does not cost more than a. At this point, we remove the 
nodes V — {W U {r}), and merge the capacity of oppositely directed edges to 
get an undirected capacity reservation x' . By construction, the pair {x' ,z) is a 
feasible solution to LPkMSriw, W, r, k) of cost at most 2 • J2e w{e)xe- □ 



We are now ready to describe our algorithm for rooted Union /c-MST. In a 
preliminary step we guess the largest distance L in the optimal solution between 
any connected node and the corresponding root, and discard nodes at distance 
larger than L from their root. This introduces a factor 0{nh) in the running 
time. Note that L < opt. We let V be the remaining nodes in layer i. 

Then we compute the optimal fractional solution OPT* = (a;', z*, z)i, of cost 
opt*, to the following LP relaxation LP^^mst for the problem, where variables 
x\ and z„ indicate whether edge e is included in the solution of layer i and node 

V is connected in layer z, respectively. Variable indicates whether node v is 
connected in at least one layer. 

min Y.i=i hEee-E^'(e)4 

E.=i,...,.<>^., "^veV-R- 

"^Ivev-R — ^' 

zl,xl > 0, 1 > z„ > 0, Vi € {!,..., /i},Vw e V - Rye e E. 

Given OPT*, we compute for each layer i a subset of nodes W'\ where 

V belongs to PF* iff zl = maxj=i^..._?i{z^} (breaking ties arbitrarily). We also 
define fc* :— [J^vgW' -^^J- each layer i, we consider the fc-MST instance on 
nodes VF'U{r*} with target fc\ This instance is solved using the 2-approximation 
algorithm apx-kmst of Lemma |8l the resulting tree is added to the solution 
for layer i. Let k' be the number of connected nodes. If k' < k, the algorithm 
connects k — k' extra nodes, chosen greedily, to the corresponding root in order 
to reach the global target k. 

Theorem 11. There is a 0{h) -approximation algorithm for rooted UNION k- 
MST. The running time of the algorithm is 0{{nh)^^^'>). 

Proof. Consider the algorithm above. The claim on the running time is trivial. 
By construction, the solution computed is feasible (i.e., it connects k nodes). It 
remains to consider the approximation factor. 

For each v G VF*, we let = Zy, and set zl — for the remaining nodes. 
Furthermore, we let x\ = h ■ x\. Observe that (i%z*,z)i is a feasible fractional 
solution to LPukMST of cost h ■ opt*. Observe also that (x%z*) is a feasible 
solution to LPksriw^, W^,V^, r^,k^): let apx^ be the associated cost. By Lemma 
[9l there is a fractional solution to LPkMSriw^ , r^,k^) of cost at most 2apx^. 
It follows from Lemma [8] that the solution computed by apx-kmst on layer i 
costs at most Aapx^ + L. 

Since the W^"s are disjoint, the algorithm initially connects at least ^ 
k — h nodes. Hence the cost of the final augmentation phase is at most h ■ L < 
h ■ opt. Putting everything together, the cost of the solution returned by the 
algorithm is at most: 

^(4 • apx' + L) + h- L<Ah- opt* + 2h ■ L < Qh ■ opt . □ 



The constant multiplying h in the approximation factor can be reduced with 
a more technical analysis, at the cost of a higher running time. We also observe 
that the integrality gap of LPuumst is Q{h) (see Appendix IB.3[) . 

3.3 Unrooted Union fc-MST 

Theorem 12. Unrooted Union fc-MST is not approximable in polynomial time 
for an arbitrary number h of layers unless P = NP . 

Proof. We give a reduction from SAT: given a CNF boolean formula on m! 
clauses and n' variables, determine whether it is satisfiable or not. For each 
variable z, we create two nodes ti and fi. Intuitively, these nodes represent the 
fact that i is true or false, respectively. Furthermore, we have a node for each 
clause. Hence the overall number of nodes is n = 2n' + to'. We create a separate 
layer for each variable i (i.e., h = n'). In layer i, we connect with an edge of 
cost zero ti (resp., fi) to all the clauses which are satisfied by setting i to true 
(resp., to false jj. The target value is fc = n' + m' . Note that, there is a satisfying 
assignment to the SAT instance iff there is a solution of cost zero to the Union 
fc-MST instance. □ 

For h = 0(1), the rooted and unrooted versions of the problem are equivalent 
approximation- wise. In fact, one obtains an approximation-preserving reduction 
from the unrooted to the rooted case by guessing one node in the optimal 
solution per layer: this introduces a polynomial factor 0{n^) in the running 
time. We remark that an exponential dependence on h of the running time is 
unavoidable in the unrooted case, due to Theorem [121 An opposite reduction is 
obtained by appending n dummy nodes to each root (distinct nodes for distinct 
layers), with edges of cost zero, and setting the target to fc -|- hn. The following 
result follows. 

Corollary 13. Unrooted Union fc-MST is APX-hard for any h> 1. There is a 
0{h) -approximation algorithm for the problem of running time 0{{hn)^^^K^). 

4 Conclusions and Open Problems 

In this paper, we introduced multi-layer covering problems, a new framework 
that can be used to describe a wide spectrum of yet unstudied problems. We 
addressed two natural ways of combining the layers: intersection and union. We 
gave multi-layer approximation algorithms, as well as hardness results, for a few 
classic covering problems (and their partial covering versions). There are several 
research questions that merit further study. 

— There are other natural ways one can combine the layers. Consider, for exam- 
ple, the car/bike problem in the case where you can put your bike in the car 
trunk. Now you can make more than one tour by bike, the only requirement 
being that the bike tours all touch the (unique) car tour. 

® Without loss of generality, we can assume that each clause does not contain both a 
literal and its negation. 



— What about min-max multi-layer problems, where the goal is minimizing 
the maximum cost over the layers? 

— We considered covering problems: what about packing problems? 

— Our algorithms for union problems give tight bounds only with respect to 
the corresponding natural LPs. This leaves room for improvement. 

— There is a considerable gap between upper and lower bounds for intersection 
problems. In particular, our hardness results do not depend on h, while the 
approximation ratios deteriorate rather rapidly for increasing h. 
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A Intersection Problems 



In this section we give the omitted details about the intersection problems. 
A.l Intersection fc-MST 

In this section we prove a simple approximation algorithm for Intersection 
/c-MST. Recall that here we are given a graph G = (V, E) on n nodes, and h 
edge- weight functions w^, . . . ,w^. By taking the metric closures of we may 
assume that G is complete. The goal is computing a tree T* for each layer such 
that Y^. w\T^) is minimized and | fl^ V{T^)\ > k. 

The algorithm is very simple: We consider a new metric w defined as a sum 
w{e) := w*(e) for each e € E, and compute a 2-approximate solution of the 
resulting (one-layer) fc-MST problem using the algorithm in 18 . 

Lemma 14. Let K QV , and w^{K) denote the cost of the minimum spanning 
tree of K on layer i. Then there exist two nodes u,v d K such that w^{u,v) < 
4w*{K)/\K\^/^ fori^l,...,h. 

Proof. Let us prove the following claim by induction on i: for any i e {0, . . . , /i — 
1}, there exist a nodeset Ki C K and paths Pi, Pi ■ ■ ■ ,Pi on Ki such that: 
(a) \K,\ > and (b) u;^(P/) < 2w^ {K)/\K\^/'' for j = l,...,i. Triv- 

ially Kq = K satisfies the claim, hence assume i > 0. Let be the min- 
imum spanning tree of K on layer i. Duplicate its edges, compute an Euler 
tour, and shortcut duplicated nodes. Let C" be the resulting cycle on K of 
length at most 2w^{K). Remove up to edges from C* so as to obtain 

jii'l"'^/'' segments of length at most 2w^{K) /\K\^/^ each. Let P be the segment 
maximizing the cardinality of Ki :— V{P) n Ki^i. Set Ki satisfies (a) since 
\Ki\ > > The paths P^ and P/, j < z, satisfy- 

ing (b) are obtained from P and P/_i, respectively, by shortcutting the nodes 
not in Ki. 

Similarly as above, we can split G'^ into \K\^^'^/2 segments which span K 
and have length at most 4:w'^{K)/\K\^^'^ each. At least one of these segments 
contains 2\Kh~i\/\K\^^'^ > 2 nodes of Kh-i. Thus there are two nodes u and v 
such that w'{u,v) < Aw'{K)/\K\^/'' for i = 1, . . . , /i. □ 

Theorem 15 (Theorem [3] restated) . The Intersection fc-MST algorithm 
above is l&k^^^^^ -approximate. 

Proof. Consider the following process: starting with the optimal set Kq of fc 
covered nodes, we iteratively take the edge {z, y} guaranteed by Lemma [T4l and 
contract it in all layers, until Ko collapses into a single node. The contracted 
edges form a tree T' (same for all layers) spanning fc nodes, of cost 

h fc— 1 h 

w{T') < 'i^w\Ko)Y^{k-i + lYT^ <A^w\Ko) / x'T^dx 

i=l 1=1 i=l 

h 

< 8k^-i^w''{Ko) = Sk^-T^opt. 

i=l 



The algorithm returns a sohition of cost at most 2w{T'). The clami fohows. □ 



We show that the bound proven m Section [52] is tight up to a factor 0{h). 

Lemma 16. The approximation ratio of the INTERSECTION fc-MST algorithm 
in Section\MMis Q[\n^~^/^) 

Proof. Take an arbitrary integer N > 2 and set n = 2'*^ — 1. We are going to 
construct weights , w'^, . . . , w'^ on an n-node complete graph G — (V, E) such 
that w*(T') = n — 1, but w{T) ~ ]7(ri^^^/''). Here T* is the minimum spanning 
tree on layer i. 

We take V — {0, 1, . . . , n — 1}, that is, the set of nodes of G are all numbers 
with up to hN digits in binary, except for 2''^ — 1, i.e., the number with hN 
ones in binary. Given x G V^, we split its /lA^-digit binary representation into h 
segments of length N and denote the i-th segment by Xi. In other words, the 
binary representation of x is x\X2 ■ ■ ■ Xh and each A^-digit binary string. 

By x{i) we denote number represented in binary that 
is, we rotate cyclically {i — l)N bits. 

To construct metric w^, sort V according to numbers x{i) and connect V 
into Hamiltonian cycle Ci in this order. All edges on Ci have weight 1 and other 
distances are minimum length distances on C^. Clearly, w'^{T^) — n~ 1. 

It is sufficient to show that, for each edge a;, y, we have w{x^ y) — Q{n^~'^/^). 
This leads to the claimed bound on w{T). We distinguish a few subcases. 

Case 1. There exists i, I < i < h, such that \xi — yi\ > 2 and {xi,yi} ^ 
{0, 2^ — 1}. Then in Gi the distance between x and y is at least n^^^/'*. 

Case 2. There exists i, 1 < i < h, such that Xi — yi. Take j such that 
Xj 7^ j/j but Xj+i = yj+i (with Xh+i — xi). Then in Gj the distance between x 
and y is at least n^^^^'^ — n^^^/'*. 

Case 3. There exists i, 1 < i < h, such that \xi — yi\ = 1. Then in Gi^i 
(with Go = Gh) the distance between x and y is at least n^^^/'^ — 2n^^^/''. 

Case 4. As 2''^— 1 is not in V, there exists i, 1 < i < h, such that Xi = 2^—1, 
yi — 0, but Xi+i = and y^+i =2^ — 1. Then the distance between x and y in 

is at least n^-^^'' ~ 2ni-2/'». 

As we exhausted all possibilities, the bound on w{T) is proven. □ 

A. 2 Approximation hardness 

We start with the following technical lemma. 

Lemma 17. Assume we have an undirected graph G — {V,E) and an induced 
subgraph G[X], X (- V with x nodes and y edges. Let 2 < xq < x and let 
c := x/xq. Then one can in polynomial time find an induced subgraph G[Y] on 
Xq nodes with at least y/(2c^) edges. 

Proof. By the linearity of expectations a random subset Xq C X containing Xq 
vertices induces a subgraph G[Xo] with — ~ 2c'' ^dges. We can 

derandomize this procedure with standard techniques. □ 



Figure 2 Construction from Lemma [T51 applied to a 4-cycle. 



Now we reduce the domain of /c-Densest Subgraph to bipartite graphs 
(see also Figure [5]). 

Lemma 18. Assume there exists a f{n, k) -approximation algorithm for fc-DENSEST 
Subgraph on bipartite graphs. Then there exists a 8/(2n, 2k) -approximation al- 
gorithm for the same problem on arbitrary graphs. 

Proof. Assume we have an instance of fc-DENSEST Subgraph, i.e., a graph G = 
(V, E) and one integer k. Construct a bipartite graph G' = {V1UV2, E') as follows: 
for each v ^ V we take two copies vi e Vi and V2 € V2. For each {u, v) ^ E we 
add {ui,V2) and {u2,vi) to E' . Let us run the /(n, fc)-approximation algorithm 
for the bipartite graph G" and k' := 2k. This way we obtain a set X' C ViU V2 
such that G[X'] has y' edges. Take X -.^ {v: vi e X' or V2 G X'}, k<\X\< 2k. 
The graph G[X] has at least y' /2 edges. Reduce X to size k using Lemma [T71 
obtaining a solution with at least y' /16 edges. 

Let us now bound how much the obtained solution is worse than the optimal 
solution. Let Xopt be any optimal solution in G such that G[Xopt\ has yopt edges 
and k nodes. In G' set X'^^^ = {wi,W2 : v £ Xopt} has 2k nodes and 2yopt edges, 
thus y' > 2yopt/ f{2n,2k). Therefore G[X] has at least yopt/{8f{2n,2k)) edges. 

□ 

Now we relate fc-DENSEST SUBGRAPH and Minimum £-Edge Coverage. A 
lemma similar to the following one was proved in |26j . but we include the proof 
for the sake of completeness. 

Lemma 19. Assume there exists a f{n)- approximation algorithm for MINI- 
MUM £-Edge Coverage on bipartite graphs. Then there exists a 2(/(n))^- 
approximation algorithm for fc-DENSEST SUBGRAPH on bipartite graphs. 

Proof Assume we have an instance (G = {V,E),k) of fc-DENSEST Subgraph. 
For the graph G we run the approximation algorithm for Minimum ^-Edge 

Coverage with consecutive £ := 1,2,..., obtaining solutions Xi,X2, We 

stop when G[Xi^i\ has more that f{n)k nodes. Assume that Xi was the last 
solution with at most f{n)k nodes. We reduce Xi to size fc using Lemma fTTl and 
return the reduced set. 




Let us now prove that it is in faet a 2(/(ri))^-approxiniation. Let Xopt be the 
optimal solution for the fc-DENSEST Subgraph instance with yopt edges. Note 
that £ > Tjopt, as was the last solution with at most f{n,)k nodes and our 
algorithm is a /(n)-approximation. Thus, by Lemma 1171 the returned solution 
has at least £/{2{\Xi\/kf) > yopt/{2{f{n))^) edges. □ 

Pipelining Lemmas [19] and [18] proves the following theorem. 

Theorem 20. If there exists an f {n)- approximation algorithm for Minimum £- 
Edge Coverage on bipartite graphs, then there exists a 16{ f (2n))^ -approximation 
algorithm for fc-DENSEST SUBGRAPH on arbitrary graphs. 

We conclude this section by showing the missing reductions for Intersec- 
tion fc-SET Cover, Intersection fc-MST and Intersection fc-METRic Fa- 
cility Location. 

Lemma 21. // there exists an f(n,k)- approximation algorithm for unweighted 
Intersection fc-SET Cover on two layers, then there exists an f{m,£)- approximation 
algorithm for Minimum £-Edge Coverage on bipartite graphs. 

Proof. Let {G,£), G = (V"i U V2,E), be the considered instance of Minimum 
£-Edge Coverage. For each u g V1UV2, let S{v) be the set of edges incident 
to V. Consider the 2-layer Intersection fc-SET Cover instance {U,k,S^,S^) 
with: U — E, k = £, and 5* = {S{v) : v e Vi} for i = 1, 2. Note that a solution 
for {G,£) translates to a solution for {U, k,S^,S^) and vice versa. □ 

Lemma 22. // there exists an f{n,k)- approximation algorithm for INTERSEC- 
TION fc-METRiC Facility Location on two layers, then there exists an f{m, £)- 
approximation algorithm for Minimum ^-Edge Coverage on bipartite graphs. 

Proof. Let {G,i), G — (Vi U V2,E), be the considered instance of Minimum 
£-Edge Coverage. For each v eViUV2, let 5{v) be the set of edges incident to 
V. Consider the 2-layer Intersection fc-METRic Facility Location instance 
defined as follows. Let C = E, k — £ and J"' = Vi. We define all opening costs to 
be equal to 1 and all connection costs w^{e, v) to be equal to if w is an endpoint 
of e, or 00 otherwise. As each client (edge) e is connected by a finite distance 
to only one facility in each layer, costs are metric. Note that a solution for 
{G,£) translates to a solution for (C, k,T^,J^'^) and vice versa. □ 

Lemma 23. // there exists an f{n,k)- approximation algorithm for INTERSEC- 
TION fc-MST on two layers, then there exists an f(n + m+ \,£)- approximation 
algorithm for Minimum ^-Edge Coverage on bipartite graphs. 

Proof. Let (G, ^), G = (Vi U 1/^2, -E), be a bipartite instance of Minimum £-Edge 
Coverage. We show how to construct the first layer of the corresponding In- 
tersection fc-MST instance: the construction of the second layer is symmetric. 

Consider the following auxiliary weighted graph G' = {V ,E'). The nodeset 
V is given by {r} U Vi U U i?, where r is a newly created root node. Moreover, 



E' = E'^U E'l^U E'^, where: E'^ = {{r,v} : v £ V2}, E'^ = {{r,v} : v G V^i}, and 
E'^ = {{v,{v,u}} : V G Vi,{v,u} G E}. We set the weight of edges in E'^^, E'^^, 
and i?^ to 00, 1 and 0, respectively. Intuitively, we want all nodes of G" that 
correspond to V2 to be too expensive to be used in this layer. Eventually, we 
consider the metric closure of G" (set E' induces a tree). We set the target k := £. 

Because of the 00 edges, no node v G Vi U V2 will belong to the intersec- 
tion. Thus the only nodes in the intersection will be the nodes corresponding to 
edges of the original bipartite graph. Consequently a solution for Intersection 
fc-MST of cost a in G" translates to a solution for the Minimum ^-Edge Cov- 
erage instance with a vertices and vice versa. □ 

Pipelining each of the Lemmas [5T1 [22] and [231 with Theorem [20) we prove Theo- 
rems [4] and [5l 

A. 3 Intersection fc-NoNMETRic Facility Location 

In this section we give a (4/c^~^/''^og^/''(/c))-approximation algorithm for Inter- 
section fc-NoNMETRic Facility Location. The algorithm works in N^'^^'i 
time, i.e., a polynomial time for any fixed h. 

Let us define the problem in a way which is convenient for our purposes. 
For each layer i, we are given a collection of subsets of the set C of clients, 
with one subset X for each facility f{X). Intuitively, X is the set of clients 
that facility f{X) is allowed to serve on layer i. We use o^{X) and w'^{x,X) as 
shortcuts for o'(/(X)) and w^{x, f{X)), respectively. A solution consists of a 
subset K C C oi k clients, and a subset C J^*. Intuitively, K is the subset 
of clients that we decide to connect, and ^xeA^ f{^) is the set of facilities that 
we open on layer i to connect K on that layer. We denote as (x) the facility 
serving client x on layer i. We slightly generalize Intersection fc-NoNMETRic 
Facility Location by adding for each x e C a cost o{x) of using x in the 
solution. In total, the cost of a solution is: 

xeK i=l XeA^ xeK i=l 

As we shall see, this generalization does not make the problem much harder, but 
allows us to describe our algorithm more neatly. 

A 0(logfc)-approximation for our generalized Intersection /c-Nonmetric 
Facility Location problem for /i = 1 follows easily from [40] . 

Lemma 24. There exists a {1 + Ink) -approximation algorithm for the gen- 
eralised fc-NONMETRiC Facility Location problem, where connecting client 
X Cz C has extra cost o(x) . 

Proof. We can reduce tha generalised version of the fc-NoNMETRic Facility 
Location problem to the classical one by simply increasing all distances between 
client X and any facility / by o{x). The claim follows from [40]. □ 



Figure 3 Approximation algorithm for Intersection /c-Nonmetric Facility 
Location. 



1: procedure FLl(fc,C, h, (J^')f^i, (o, w,o)) 

2: it h^l then 

3: return solution found by Lemma [24] 

4: if fc = 1 then 

5: return optimal solution by brute-force 

6: if ^ 0, A' for i = 1, 2, . . . , h. 

7: repeat 

8: for all 1 < r < ft and all X G J^'' do 

9: for j — 1 to min(fc - lA"!, IXVi^l) do 

10: o' -s— o, except for elements x £ X \ K , where we put o'{x) = o{x) + 
w*(x,X). 

11: iKj,x, iAlx)a=i) ^ FLIO', X\K,h~l, {F%<,<H,i_^r, {o,w,o')). 

12: Cj,x cost of {Kj^x, {A)^x)^=i) w.r.t. costs (o, w,o'), plus o'^{X) 

13: c,,x ^ Cj,x/i. 

14: ro,Xo, jo values of the loops' iterators for which the cheapest solution 

was found according to weights Cj,x- 

15: K^KU Kj„,Xo , A"-" ^ A"-" U {Xo}, A' ^ A' U A^^^Xo ^r i + ro 

16: For each x £ Kj^^Xg assign A'^° (x) = Xq. 

17: until 17^1 = k 

18: return {K, {A')ti) 



Now we solve the problem for arbitrary h. Let us give some notation. We 
are going to develop a procedure FLl(fc,C, h, (J^')jLii {o, w,o)) that returns the 
(4fc^~^/'' log^/''(fc))-approximation for Intersection fc-NoNMETRic Facility 
Location. We abandon the requirement that each X £ U!Li is contained in 
C, but the solution set K needs to be a subset of C. 

Note that Lemma [Ml provides an algorithm for /i = 1. For arbitrary h, we 
use the procedure described in Figure [3] The FLI procedure constructs set K 
with facilities and a choice function greedily. At one step, we iterate over all 
layers r and sets X £ J"^ and all possible cardinalities j of elements that set 
X can cover (i.e., 1 < j < min(fc — \K\, \X \ K\)) and try to: choose X and 
solve problem without sets J-^ (where X £ T^) for universe restricted to X \ is' 
and the goal size j. We solve the subproblem using FLI procedure, but for h 
decreased by one. We hide the costs of attaching elements x £ X \ K to X in 
costs per element, i.e., in o{x). Finally, we choose set X and cardinality j with 
the smallest cost per element covered. This cost is kept in variable Cj^x- 

Proof. (Theorem\^ We consider the algorithm in Figure [3] It is clear that FLI 
procedure works in A^'^C') time, as it calls recursively itself 0{Nk^) times with 
h decreased by one. 

We next prove the claim on the approximation by induction on h. For h = I 
the thesis is implied by Lemma [Mj Assume now that all recursive calls in the 
FLI procedure return solutions with (4fc^~^/(''~^^ log^^''''~^''(fc)) approximation 
ratio. 



Let opt be the cost of the optimal solution for the given instance. Pick any 
optimal solution with cost opt and let Ko be the set of covered elements by it, 
and be the chosen subset of J"' for i = 1, 2, . . . , h. For each layer i and each 
element x e Ko we fix the set (x) e that covers x in the optimal solution. 

We prove that at one step of the algorithm, the weight cj^^Xo satisfies: 

cj„Xo < opt ■ 4i-i/^(fc - lATI)-!/^ \og'/\k). 

Recall that Cj^^Xo is the average cost paid newly covered elements. This bound 
is sufficient, since the total cost of the constructed solution is bounded by: 

total cost < opi • V 4i-i/''(fc - i)-^/'' log'^'^ik) < opt ■ ^^'^'^ \oi'^{k) [ x-^'^dx 
= opt ■ logi/'^/t)^— i-^fci-i/'^ < opi • 4 \og^''^{k)k^-^/^. 

The last inequality follows from the fact that 4~^ < 1 — eforO<£< 5. 

Let K := A^~^/^{k— \K\Y~^/^ \o^^^{k). We consider two cases, depending on 
whether there exists a layer r and a set X e that covers at least /t elements 
of Ko \ K, i.e, 

\{xGKo\K: C {x) = X}\ > K. 

Case 1. Assume there exists a layer r and a set X E C such that for at 
least K elements x of Ko \ K we have C (x) = X. Let us focus on the mo- 
ment when our algorithm considers taking set X. We may assume k < k — \K\, 
as otherwise k — \K\ is boundcid by constant and we may instead use brute 
force to finish the greedy construction optimally. Therefore our algorithm con- 
siders covering k elements of X. As the optimal solution does it, it may be 
done with cost opt, so the recursive call returns the solution with cost at most 
opt ■ 4k;^~^/('*~^) log^^'-''~^^(K). We cover k elements, so 

Ck,x < opt ■ 4k-V('»-i) log^/e^-^^C^) 

-1 

<opt-4 (4^(fc - \K\)^ logi/^fc)) logi/e^-i^fc) 
= opt ■ 4i-i/''(fc - \K\)-'/''\og^/''{k). 

Case 2. Every X G T'^^^ covers at most k elements of Ko \ K. For each 
X e Ko \ K denote 



h 

W{x) = o{x) + w''{x, O'' (x)) + o^'iO'' (x)), 



i.e., the total cost of choosing x, attaching it to set C" (x) and choosing set 
{x). By the assumption in this case, we have 

h 

xeKo\K xeKo\K r=l 

h 

+ ^''(^) -lix & Ko\K -.O'' {x) =X}\<K- opt. 

r=i xec- 

Thus there exists xo & Kq \ K such that w{xq) < k ■ opt/\Ko \K\. Let Xq = 
(xo). Note that, since our algorithm uses brute force for fc = 1, the recursive 
call with find optimal solution for j := 1 and X := Xq and thus Ci^Xo ^ w{xo). 
As\Ko\K\>k-\K\we\i&Ye: 

c\,Xo < K • opt/\Ko \ K\ 

< opt ■ 4i-i/''(fc - logi/''(fc)(fc - \K\)-^ 

= opt ■ 4i-i/'*(fc - \K\)-'^'' logi/''(fc). 

□ 

B Union Problems 

In this section we give the omitted details concerning union covering problems. 
B.l A Greedy Approach 

We next describe a simple greedy algorithm which provides a logarithmic approx- 
imation for several union partial covering problems. Consider a partial covering 
problem whore U is the set of requests, 5* is the set of items on layer i, with 
costs : <S* — >■ R>o, and k is the target. We require that the covering problem 
satisfies a natural composition property, namely two solutions satisfying k' and 
k" distinct requests, can be merged (without increasing the total cost) to obtain 
a solution satisfying k' + k" requests. (Merging might involve some polynomial- 
time operations). The algorithm works as follows: 

(1) For all layers i, for all fc* := 1, . . . , fc, solve the single-layer problem induces 
by the triple {U,S^, fc') with a ^-approximation algorithm. 

(2) Among all the solutions computed, take the one A, obtained for some triple 
{iU,S'^,k'^), which minimizes the ratio w^{A)/k'^. 

(3) Merge A with the solution under construction. Remove from U the requests 
satisfied by A, and decrease fc by fc*. 

(4) If fc > 0, go to Step (1). Otherwise return the current solution. 

Theorem 25. The algorithm above computes a 0{p log k)- approximation for 
the partial covering problem considered in polynomial time. 



Proof. The claim on the running time is trivial. The algorithm computes a feasi- 
ble solution, due to the composition property. Consider now the approximation 
ratio. Let Ai, . . . ,Aq be the sequence of approximate solutions computed, wj 
be the cost of Aj and kj the number of requests that it satisfies on layer ij. 
Observe that at the beginning of iteration j, the current number of requests is 
k — X)o<j and the cost of the optimal solution with respect to that number 
of requests is no more than opt. By an averaging argument, at each iteration j 
we have w^^{Aj)/kj < opt/{k — J2a<j ^i)- '^^'^ conclude that the cost of the 
solution computed is at most 



Corollary 26. There are 0(log k)- approximation algorithms for Union /c-MST 
and Union fc-METRic Facility Location. 

Proof. Observe that removing requests transforms the original fc-MST problem 
in each layer into a fc-STEiNER Tree problem: for the latter problem there is a 
4-approximation algorithm [18 . Note also that all the partial solutions in layer 
i contain the root r*: hence the merging step is trivial. The claim for Union 
fc-MST follows. 

For /c-Metric Facility Location, there is a 2-approximation algorithm in 
|29| . In this case removing a request simply means removing one client, and the 
merging step is trivial. This proves the claim for Union /c-Metric Facility 
Location. □ 



B.2 Union fc-METRic Facility Location 

In this section we present an LP-based 0(/i)-approximation algorithm for Union 
fc-METRic Facility Location. As we wiU see, the basic idea is an for Union 
fc-MST. 

RecaU that in Union fc-METRic Facility Location we are given a graph 
G = {V,E), a set C C U of clients, a set C V oi facilities, one integer fc 
(target), a set of opening cost functions o^ : T ]R>o, and a set of edge-weight 
functions : E ^ II^>oj with i = 1, . . . ,h. The distance between nodes u and v 
w.r.t. is denoted as w'^{u, v). A feasible solution is given by a pair (C*,.F') for 
each layer ?, C C C and J^* C such that | Ui C'| > fc. The goal is minimizing 
the cost X]i=i hi^fef^ -I- X^csc* -^*))- Here w^{c,J^') denotes the 
minimum distance on layer i between client c G C and facility f ^ T' J- . 




□ 



Also in this case we consider a natural LP relaxation LPkMFL{C, o, w, k) 
for the single-layer version of the problem: 



min E/e^o(/)2// +X;(c,/)ecx^^(c,/)a:c,/ 

si- Xcj < Vf, V(c, f) eC X J^; 

I^/e^2;cj > Zc, VceC; 

a;c,/,y/ > 0,1 > > 0, VceC,V/GJ-. 

Variable indicates whether facility / is opened, and variable Xcj whether 
client c is connected to facility /. Variable Zc indicates whether client c is con- 
nected to some facility. We need the following result. 

Lemma 27. |77f Let {C,T,o,w,k) be an instance o/fc-METRiC Facility Lo- 
cation, Omaa; := maxygjr o(/), andopt' be the optimal solution to LPkMFLiC,J'TO,w,k). 
There is a polynomial time algorithm apx-kmf 1 which computes a solution to the 
instance of cost at most iopt' -f- 2omax ■ 

The algorithm that we use is analogous to the one for the fc-MST case. In a pre- 
liminary phase we guess the largest cost o* of a facility in the optimum solution, 
and remove all facilities of larger cost. Let be the remaining set of facilities 
on layer i. We then compute the optimal solution OPT* = (x*, z*, z)^, of cost 
opt* , to the following relaxation LPuumfl for the problem: 

min E»=i,...,/.(E/G^. o'{f)y} + E(c,/)GCx^» /) </) 

s.t. xlf<y), Vie{l,...,/i},V(c,/)eCx J-'- 

E/g^.</>^c, Vje{l,...,/i},VceC 

E.=i,...,/.4>^c, VceC 

a;^ J, y}, > 0, 1 > > 0, Vf e {1, . . . , /i}, Vc e C, V/ e ^. 

Then we identify for each layer i the subset of clients C* := {c e C : z^ = 
maxjjz^}}. We run the algorithm apx-kmf 1 from Lemma [27] on each layer, 
with clients C*, facilities and target fc* := EcsC' •^^J: we open facilities and 
connect clients accordingly. Let k' be the number of connected clients. If k' < k, 
we connect extra clients in a greedy fashion, possibly opening new facilities: in 
particular, we consider the pairs (c, /) € x T^, with c not connected, which 
minimize o*(/) -|- w'{c, /), and we connect the corresponding clients. 

Theorem 28. There is a 0(h) -approximation algorithm for UNION fc-METRiC 
Facility Location. The running time of the algorithm is 0{{nh)^^'^'>) . 

Proof. Consider the algorithm above. The claim on the running time is trivial. 
As in the fc-MST case, consider the feasible fractional solution {x^,y^, z% z)i ob- 
tained from OPT* by setting z^ = Zc if c S C*, 5J = otherwise, and raising the 



variables x and y by a factor h. This new solution costs at most h ■ opt* . Further- 
more, {x^,y^,z^) is a feasible fractional solution to LPkMFL{C^,J'^,o\w^,k'^). 
Let apx^ be its cost. Lemma \T7\ guarantees that the cost of the integral solution 
on layer i is at most 3apx^ + 2o* . Since J^i k'^ > k — h, the final step costs at 
most h ■ opt. Altogether the cost of the solution computed is at most 

^(3apa;* + 2o*) + h ■ opt < 3h ■ opt* + 2h ■ o* + h ■ opt < 6h ■ opt. 

i 

a 

Also in this case a more technical analysis allows one to reduce the constant in 
front of h in the approximation factor, at the cost of a larger running time. 

B.3 Integrality Gap 

Lemma 29. The integrality gap of LPukMST cind LPukMFL is f2{h). 

Proof. We consider the following unweighted Set Cover instance given in [44] . 
Let G' be an hypergraph on m! nodes, which has one hyperedge for any subset 
of m'/2 nodes. We construct a set cover instance with m' sets given by nodes, 
(m72) elements given by hyperdeges, and inclusion given by incidence. Taking 
a fraction 2/m' of each set gives a feasible fractional solution of cost 2 to the 
natural set cover LP. On the other hand, the optimal integral solution uses 
m'/2 + 1 sets. Hence the integrality gap in this case is fi{m'). 

The same reductions as in Theorem [7] imply a fi{h) lower bound on the 
integrality gap of LP^kMST and LPukMFL for the case k = n. □ 



